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Abstract 

O 
in 

I Process of quantum tunneling of particles in various physical systems can be 



effectively controlled even by a weak and slow varying in time electromagnetic 
signal if to adapt specially its shape to a particular system. During an under- 

e , 

barrier motion of a particle such signal provides a "coherent" assistance of 

£3 . 

tunneling by the multi-quanta absorption resulting in a strong enhancement 
of the tunneling probability. The semiclassical approach based on trajecto- 

H : 

ries in the complex time is developed for tunneling in a non-stationary field. 
Enhancement of tunneling occurs when a singularity of the signal coincides 
in position at the complex time plane with a singularity of the classical New- 
tonian trajectory of the particle. The developed theory is also applicable 
to the over-barrier reflection of particles and to reflection of classical waves 
(electromagnetic, hydrodynamic, etc.) from a spatially-smooth medium. 
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I. INTRODUCTION 



Control of quantum systems by tuned external signals is an actively developed field 
at present, see for example |IJ and references therein. Excitation of molecules, when one 
should excite only particular chemical bonds 0-0, formation of programmable atomic wave 
packets [H, a control of electron states in heterostructures 0, and a control of photo-current 
in semiconductors J7| are typical examples of control by laser pulses. A control of quantum 
tunneling by electromagnetic signals is also a matter of interest, since tunneling is a part of 
many physical processes and of some chemical reactions. Modern facilities enable to tailor 
fast signals required for this purpose |||| . 

Let us focus on main aspects of tunneling under a non-stationary perturbation. The po- 
tential barrier V(x), extended over the distance a, sets two typical energy scales: the barrier 
hight Vq and tkv ~ h^Vo/ma 2 , where u can be associated with an oscillation frequency in 
the overturned potential. For semiclassical barriers the two energy scales are well separated 
hu <C Vq and in absence of a non-stationary field the probability of tunneling through the 
barrier can be estimated with the exponential accuracy as W ~ exp(— Vo/hu). In presence 
of the periodic signal Sq cos fit a particle can absorb the quantum Q with the probabil- 
ity (aSn/Vo) 2 and tunnel in the more transparent part of the barrier with the probability 
exp(— (Vo — Q)/hw). The total tunneling rate can be written as 

— B)+(t)-(-") 

Eq. 1 is approximate since in quantum mechanics one should multiply amplitudes but not 
probabilities, nevertheless the form (1) accounts necessary physical feature. Suppose a 
tunneling particle is acted by some electric field £(t) (periodic or pulse- type) and the typical 
time scale of this signal is 9. Than one can distinguish two different physical situations: 
(i) hard signal, fast varying field with 9 ~ H/Vq and (ii) soft signal, slow varying field with 
9 ~ uj^ 1 . Tunneling can be governed easily by a hard signal even when its amplitude is 
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less than the static barrier field Vo/a, since the probability of quantum absorption (aS/Vo) 2 
competes, according to Eq. 1, with the small tunneling rate (in this case Q ~ O^ 1 ~ Vq/H). 
It is also obvious that a soft signal of the high amplitude Vo/a (static field of the barrier) is 
able to govern tunneling. Can a soft signal with the amplitude much smaller than the static 
field of a potential barrier control effectively a tunneling process? Suppose a soft signal has 
the following shape 

«<«> = jrT¥j¥r (2) 

Its Fourier harmonics Sq ~ (Q9) n ~ 1 £6 exp(— Q9) should be inserted into Eq. 1. As follows 
from Eq. 1, when the signal width 9 is less than l/2u, the quadratic in S correction to 
the static probability diverges with increase of Q. It means the perturbation theory with 
respect to a weak non-stationary signal to break down for sufficiently short pulses. Note, 
the pulse still remains soft. This is an indication of efficiency of soft signals. Breaking of 
the perturbation theory means a significance of multi-quantum processes and a principal 
question is that what theory has to be used in this case. 

A review of some aspects of tunneling in complex systems, including the instanton ap- 
proach, was done in the book [ nj , see also Ref . |]TTJ . Recent achievements in the semiclassical 



theory under stationary conditions are presented in Refs. ||12|-|14jl. As it has been argued 



in Refs. [|T|,|]l|[l(]] , the semiclassical method of complex trajectories is applicable also to 



a non- stationary case, when a signal is periodic in time. Nevertheless, despite a number 
of publications, use of semiclassical theory for tunneling in a non-stationary field remains 
non-obvious. What happens in general case, in particular, for a short pulse like one given 
by the relation (|2|) ? The goal of this work is to show that the semiclassical theory based on 
the concept of the complex time is an appropriate description of tunneling under action of 
a soft pulse of any shape. For the particular case of a triangular barrier the tunneling rate 
in presence of a non-stationary (soft) field is found to be determined, in the main exponen- 
tial approximation, by the classical action S(x, t) satisfying the Hamilton- Jacobi equation. 
In this approximation the wave function is proportional to exp(iS(x,t)) (Below Planck's 



constant is unity). The first correction (preexponential) to this classical result and the sec- 
ond one were found explicitly and shown to be small comparing to the main contribution, 
what is typical for semiclassical approximation. The tunneling rate is found as a function 
of time, it tends to its static value at t — > ±00, when £(t) = 0, and reaches the maximum 
at some moment of time. This maximum value is given, with an exponential accuracy, by 
an extreme value of the classical action, which is determined, according to classical me- 
chanics, by means of classical trajectories. The classical trajectory obeys Newton's equation 
md 2 x/dt 2 + V'(x) = £(t) in the complex time, since in the real time there is no classical 
under-barrier path. 

So, the method of classical trajectories in the complex time can be used, when the full 
time dependence of the tunneling rate is not required, but only its maximum value with 
an exponential accuracy is a matter of interest. Note, despite the classical action depends 
functionally on trajectories defined in the complex time plane, the argument t of the action 
S(x, t) is considered to be always real. The complex time has no physical meaning, it is only 
a possible way to parameterize a solution of the Hamilton- Jacobi equation. Nevertheless, 
this way is very useful. For a weak non-stationary field S(t) the classical trajectory Xo(t) 
satisfies the unperturbed equation (m/2)(dx /dt) 2 + V(x ) = E, where E is the particle 
energy, which can identified with an incident energy of particle flux on the barrier V(x). 
The classical trajectory Xo(t), as a function of the complex t, has the singularity at t — t s (E) 
l^jTBH and the external signal (0) has the singularity at t = %Q. As shown in this paper, 



when the two singularities coincide 

lmt s (E)=e (3) 

the effect of an external signal on tunneling enhances. Under the condition (|3]) the perturba- 
tion theory breaks down at essentially weaker non-stationary signal comparing to a general 
case. The parameter lmt s (E) depends on properties of the static barrier and the particle 
energy [jD^jnj, but 9 is a characteristic of a non-stationary signal. The physical meaning 



of the condition (H) is not straightforward since a quantum mechanical process is described 



by a product of amplitudes but not of probabilities. The condition (^|) rather corresponds 
to a coherent cooperation of tunneling and quanta absorption, in other words, it is some 
"resonance" condition between motion of system and external signal. Eq. |3| determines some 
remarkable threshold energy Ex- As shown below, when the particle energy is big E > Ex 
the tunneling process is moderately violated by the signal (0) (of course, if the signal is less, 
than the static field of the barrier V/a); for lower energy E < Et the process of the barrier 
transition is strongly stimulated even by a relatively small signal. This opens a possibility 
to manipulate effectively a tunneling process by a specially adapted (according to Eq. [D 
electromagnetic signal of a small amplitude. 

In Sections II- VII the tunneling probability as a function of time is calculated for the 
triangular barrier. In Section VIII the method of complex trajectories is described. In 
Sections IX-X this method is applied to the triangular barrier to compare with the results 
obtained by the direct solution of the Hamilton- Jacobi equation. In Sections XI-XIV the 
method of complex trajectories is applied to a barrier given by an analytical function V(x), 
when there is no simple solution of the Hamilton- Jacobi equation. 



II. TRIANGULAR BARRIER 

In this Section we consider decay of the metastable state in the potential 

V(x) = V - £ I x I -\ 2{y ~ ^ S(x) (4) 

V m 

under action of a non-stationary electric field S(t). In the limit £ o —>■ the energy E 
corresponds to the bound state in the 5-function potential well. The symmetric wave function 
(ip(x,t) = ip(—x,t)) can be written down in the form 

ip(x, t) — exj)(iS(x, t) + ia(x, t)) (5) 

where the classical action S obeys the Hamilton- Jacobi equation at x > 

dS 1 fdS\ 2 TA „ 
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with the boundary condition 

At x = one can impose the condition S(0,t) = —Et. The equation for a has the form 

da 1 OS da 1 fda\ 2 i d 2 a _ i d 2 S 
dt m dx dx 2m \dx J 2m dx 2 2m dx 2 



with the boundary condition 



Equations @ and @ are exact ones. The solution of the Hamilton- Jacobi equation (|D can 



be found by conventional methods [17 



S(x, t) = - ^- jf dh (p + £ h + jf dt 2 £{t 2 )^j 

+ (p + £ t + dt^it^ x + {V- E)t - Vt (10) 

The functions p(x,t) and t (x,t) have to be defined from the conditions dS/dp = and 
(95/ <9i — what gives the following expressions 



rk) 

p(x, t) = i^/2m(V — E) — £ t - / dt 2 £ 

Jo 



(*2) (11) 



and 



mx = J dti {p + £oh + ^ dt 2 £{t 2 )^j (12) 

Eq. [12] has to be inserted into Eqs. [11] and 13, what results in the final expression for the 
action 

S{x, t) = - J- J dt x (iy/2m(V - E) + (h - t )£ + jf" dt 2 £(t 2 )^j 

+ x ^v/2m(V-E) + (t - t )^ + J dtiSih)^ +(V- E)t - Vt (13) 
where the function t (x,t) is given by the equation 

i{t - t ) ^2m{V -E) + (t- t ) 2 ^ + fdhit-tJSfa) (14) 

1 J to 
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mx 



By means of the relation (a partial derivative is taken under the fixed lower index) 

m) x + m{fa)Xfa) t = (atX (15) 

the equation (^) for a in new variables to and t has the form 

fda\ 1 (da\\ i (d 1 da\ 1 + h(t ) 
I ~n7 ) " 7777 77777 777" + 7777 7777 I 7u B aT ) = ~ T7C ( 16 ) 



\dt J t A(V — E)F 2 \dt Q J t 4(V-E)F \dt F dt J t Ft, 



'00 

Here new notations are introduced 

F(t , t) = l + i-° (1 + h(t )) ; = 4^ r oo = 7 " 

TOO <^0 <^0 

(17) 

In terms of new variables the boundary condition (^) reads 

'da(t ,t) 



m ) =0 (18) 

In semiclassical approximation o should be small comparing to a big classical action S and 
it can be expanded into a series 

o- = o-i + o- 2 + ... (19) 

which is produced by Eq. [T^, where the last two terms in the left-hand side have to be 
considered as a perturbation. Now one can write 

a n (t ,t)= [ d77$ n (77,to)+ /°dfi$n(0,ti) (20) 
«/ o Jo 



where 



and 



- to, to) = TYTTfr" < 21 > 



^ If ® a i \ 2 i f d 1 du\ \ 

$2(t ~ to ' to) = Aty - E)F* \WJ t - A(V - E)F [df Fd^) t (22) 



From Eqs. ^ and |21] one can obtain an explicit expression 



itT 1 (t Ql t) = ~]nF(t Q ,t) + -^- I ° dt x (1 + h{t x )) (23) 
2 2t 00 Jo 
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The expression for 02 can be easily obtained from Eqs. ^ and |22] but it is too bulk and we 
do not write it here. The main parametric estimate at n > 1 

i>&n ~ 7T-T (24) 

characterizes Eq. [19] as a typical semiclassical series since (V — E)t 00 ^> 1. The pulse width 



is supposed to be order of r 00 . Eqs. 13 and 111 for the classical action and Eqs. [19| - |23| for cr 



enable to consider a decay of the metastable state under action of the non- stationary field 
£{t). 

III. CAUSALITY 

Suppose a pulse of the electric field has the form 

£{t) = S(t)e(t' -t) + £'{t)Q{t - t') (25) 

The equation (|25|) can be considered in the complex time if to represent 0-functions in the 
form 

With the definition ( ^6|) the 0-function can be treated as 0(Ret). As follows from here, 
the function £(t) at the complex t-plane depends on £{t) {£'(t)) only to the left (right) of 
the vertical line Re t = t' . If to chose the contours of integration in Eq. 13 to the left of 
the vertical line with the real part t, than the action S(x, t) does not have an information 
how the non-stationary pulse behaves at moments later than t. This choice of contours of 
integration corresponds to the causality principle. 

IV. CLASSICAL ACTION 



In this Section we consider only the classical action S in Eq. [|. Under the signal (0) 
the imaginary part of S reaches its minimum value at some moment of time resulting in a 
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maximum of the decay rate. For this reason we consider first the action S(x, 0) ait — 0. For 
a pulse symmetric in time to(x, 0) is an imaginary value and we introduce t (x) = —Uq(x, 0). 
At small amplitude of the signal (0) essential values of tq are close to 9 and the new variable 

-l-f (27) 

is convenient (z <C 1). Below only integer values n > 3 in Eq. |2| and 9 < t 00 are considered. 
It follows from Eqs. 13 and |14| 

(r 00 - 9) 1 - — £ ; zi 



(28) 



and 



. <9 2 S(x,0) _ 9/(r 00 -9) + (z 2 /z) n m _ ( 89 \ 1/n 
1 o3 - 1 /.. /_w 7T i *2 - TwZ TrTF" 



«9x 2 l-(^A) n \2"(t 00 -9)£ 



In Eqs. |8| and ^9] the amplitude of the signal £ is supposed to be small leading to small z\ 



and z%. As follows from Eq. [14|, the function tq(x) is determined by the relation 

dmx £9 , , 

By means of Eqs. |28] and |29| one can find the coordinate dependence of the action S, which 
is shown in Fig. [I]. Two branches in Fig. [I] in the limit £ = go over into conventional WKB 
wave functions exp(± J | p | dx). At the lower branch, where dS/dx = 0, 



32 \ _ £ n # 2 

3r 2 J ' Xl ~ 2m 



iS(x 1 ,0) = (V-E)9[l- Tir ) ; ^i = ^zr ( 31 ) 

r oo, 



At the common point, where d 2 S/dx 2 — ► oo, 

i£(x 2 ,O) = (V-£)0(l-— ^ ; x 2 = ^(2r o-^) (32) 

V r oo / ^ 

Near the common point x 2 the deviation of the action from the value fl3"2|) is proportional 
to (x 2 — and in the second order to ±(x 2 — x) 3 / 2 , what develops the two branches. 
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V. NON-SEMICLASSICAL CORRECTIONS 



For validity of the semiclassical approximation the inequalities 

| S |>| 0"i |>| a 2 |> ... (33) 

should hold. Like in a static case, one can expect a violation of the semiclassical theory near 
the point x = x±, where dS/dx = 0, and the point x = x 2 , where d 2 S/dx 2 — > oo. Let us 
compare o with the classical action S near these "dangerous" points. 

n/(n+2) 



Under the condition zi<z<z 2 , as follows from Eqs. 21 -12 



and 



. . i 3n(n+l) + n(2n-3)(V;)" p 2 \»+ 2 



An x-dependence in the right hand sides of Eqs. |34] and |35] comes through z according to 
Eq. At the point x — x\ (z — Zi) 

1 , fn- 1 / 9 \\ 1 m 



and 



'■^■'■j.o) -^7-—— --— — (s(//-Di i-— y+ n ( 2n 3 ) 



48(^-^)^(1-^00) \ v J \ r 00 J n-1 

(37) 

A particle positioned in the well corresponds to 5(0, 0) at the upper branch in Fig. 1 {z ~ 1). 
A particle leaves the barrier under the condition dS/dx = at the point x — x\ of the lower 
branch in Fig. 1 (z — Z\). One should be sure the points x = and x = X\ relate to the same 
semiclassical solution, on the other words, it should be possible to find a way from to X\ 
with no violation of the semiclassical conditions (|33"D. Between real points z ~ 1 and z — z\ 
there is only one "dangerous" point z = z 2) where, according to Eq. B5L a — ► oo and the 



condition (|33|) breaks down. Nevertheless, the semiclassical approximation remains valid if 
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the condition (|33"D holds on some contour | z — Z2 |~ £2 around the point Z2 in the complex 
z-plane. The point Z2 (x = X2) is a branch point of the action, which has a contribution 
proportional to (22 — x) 3//2 , like a turning point in a static problem. The sequence of the 
Stokes and anti-Stokes lines |19[] going from this point is qualitatively the same as in a static 
case. The condition ( |3"3"D on the above circle can be written in the following approximate 
form if to put z ~ z 2 in Eq. |35] {6 < r 00 ) 

n \ n/2— 1 c 



In the relation ( |38D the numerical coefficient a n ~ 1 at n ~ 1, but at big n the coefficient 
a n increases, what prescribes to choose a not big n for validity of the semiclassical approxi- 
mation. The condition | 02 |<C| <J\ \ at the point x = x% is less rigorous. The semiclassical 
conditions (^) require the pulse amplitude £ to be not small. Remarkably, this amplitude 
can be still less than the static barrier field Sq. At lower S, than one satisfying the relations 
(|38|), one should expect the perturbation theory to be applicable. 

VI. FINITE TIME 

Eq. [14] determines the function £0(^5 1) an d for the signal (fj) one can write at x = X\ and 
\t\<6 

i ^- t) = iT ' + 2(^)(4 + ^) (39) 

According to the causality principle, the contour of integration in Eq. 13 should be to the 
left of the time t and hence the condition Re to < t results in the restriction t > 0. At t < 
the semiclassical approach in its present form is not valid since the integration penetrates 
"in the future" and this case requires further investigation. The time dependence of the 
classical action can be found from Eq. 10 

2 



Im = __L i m ^2m(V-E) + (t - t )£ + jf" dhS^ 



(40) 
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At x = x\ and t = the right hand side of Eq. [40] is zero. One can easily obtain 



t 4 



Im 3(x h t) = Im S(x h 0) + (n - 1)(V — E) (41) 



Or 2 



VII. THE TRANSITION PROBABILITY 

Suppose w(t) is the probability to find a particle in the 5-function potential and initially 
w equals unity. The continuity equation reads dw/dt = — (2/m)lm(ip*dip/dx), where the 
right hand side is taken at x = x\. Using expression for the wave function 

ip(x, t) = ip(0, t) exp (iS(x, t) - iS(0, t) + i<r(x, t) - ia(0, t)) (42) 

where 

V>(0, t) ~ (2m(V - E)) 1/4 exp(-iEt) (43) 

and by means of Eqs. ^2] and ^ one can obtain 
dw(t) l8(V-E) fdReS 



(^exp(-2Im(S + a))^ ^ (44) 



dt V m 

Eqs. p | and 13 give dR.eS/ dx = £$b at x — X\ and the main time dependence in the 
exponential of Eq. [|4] follows from S(xi, t) (Eq. ^l]). Collecting the all terms, one can obtain 
finally (9 < r 00 ) 

dw(t) 2{V-E)t ( 89 N 1 /^- 1 ) / (V-E)t 4 



dt e (n-l)-/(n-i)( roo _0) \ 2 (r 00 -e)S J e ^ [ 2 ( n 1) ^ 

exp (-2(y-^(l-^)) (45) 

As discussed in Section VI, Eq. 45 is valid only at t > 0. The semiclassical conditions ( |33"|) 
are supposed to hold. According to Eq. 45, the typical time scale of the output flux At ~ 
(9tq /(V — E)) 1 / 4 determines the uncertainty of energy of outgoing particles AE ~ At^ 1 
which is much smaller than the energy E. The decay rate dw(t)/dt tends to its static value 
at t — > ±oo having a maximum at some moment of time 



dw\ ( . / 9 



J2 



^-•n-^-n 1 -ss.il (46) 
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The exponent in Eq. |4(], according to its derivation, is a minimum value of the imaginary 
part of the classical action. A minimum value of action can be calculated, as known from 
classical mechanics, by means of trajectories satisfying Newton's equation. So, when we are 
not interested in the full time dependence of a decay rate but we need only its maximum 
value with an exponential accuracy, the method of classical trajectories can be used. This 
method is described in the next Section. 

VIII. METHOD OF COMPLEX TRAJECTORIES 

In this Section we consider penetration of incident particles through a potential barrier 
under action of a non-stationary pulse. We restrict ourselves only by the main exponential 
approximation when one can use the semiclassical expression for a wave function ip(x,t) ~ 
exp(iS(x, t)). We consider here a particle flux on to the barrier shown in Fig. 2, but the final 
result can be easily applied to decay of the metastable state through the triangular barrier 
The maximum value of the outgoing flux of particles can be calculated as a maximum 
with respect to time 

W max ~ max | exp(iS(x,t) -iS(x ,t)) | 2 (47) 

Here x is some coordinate to the right of the barrier, x ~^ — oo, where lmS(xo,t) = 0, and 
in Eq. [|7| one can put S(x ,t ) instead of S(x ,t). The right hand side of Eq. f|7] does not 
depend on x and xq. It is a function of t only. Eq. |47] corresponds to the extreme classical 
action, which can be found by method of classical trajectories x(t) defined in the complex 
t-plane, since in real time there is no classical trajectory for an under-barrier motion. The 
complex path C is shown in Fig. 3. The real classical turning point is x± = x(ti), where 
dx(t)/dt = 0. This point corresponds to the classical exit of a particle from under the 
barrier. The real coordinate xq = x(to) is defined under the condition Re to = to ~~ * — oo. 
The classical trajectory connects the points {xq, to} and {x\, ti}. The contour C is symmetric 
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with respect to the real axis. It is convenient to write Eq. ^7| in the form 

W max ~| exp(iS(xi,ti) -iS(x ,t Q )) | 2 (48) 

Since at t — > — oo the non- stationary field £ (t) — > 0, a connection between values of the 
action at the points t and t Q is simple 

S(x , t ) = S(x , i ) + (*o - t )E (49) 



where E is the energy of an incident particle. According to Eqs. ^8] and [J9], the maximum 
amplitude value of the outgoing flux of particles has now the form 

W max ~ exp(-A) (50) 

where 

A = -i J c dt h (^j 2 - V(x) + x£(t) + E^j (51) 

is defined by means of trajectory satisfying Newton's equation m d 2 x/dt 2 + V'(x) = £(t) in 
the complex time. Due to symmetry of the contour C the value of A is real. The trajectory 
x(t) should not be necessary real at all t, it should be real at least in vicinities of real points 
Xq and X\. On the left horizontal parts of the contour C, where £(t) = 0, x(t) satisfies the 
equation 



m f dx\ 



2 



2\at) +v{x) - E (52) 

and is expressed through the real functions x = f(t — t , E) (up) and x = fit — i*^ E) (down), 
where E is a real energy. Now one can formulate conditions how to choose the contour C: 
for given (at t — > — oo) particle energy E and the pulse shape £{t) one should find Imt 
and the real turning point X\ = x(ti). Eq. [51] holds for a potential barrier V(x), which is an 
analytical function of the variable x. Such a barrier has no artificial restriction in coordinate 
(no singularity at a real x). For this reason, the equation fl5"T| ) can be interpreted as one 
accounting not only an under-barrier part but also some pre-barrier motion. 
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IX. APPLICATION TO A TRIANGULAR BARRIER 



Eq. BT] is applicable to the case of a potential barrier V(x), which is an analytical function 



of the variable x. In the case of the triangular barrier (|j), which is a non-analytical function, 
the all classical path of the particle is restricted by an under-barrier motion. In this case 
the exponent A, instead of Eq. [51], should be written in the form 



A = 2 Im 



jf dt (? (If) 2 _ v + xSo + x£{t) + E ) (53) 

where x(t) is the trajectory satisfying the Newton equation in the complex time 

m^-£ = S(t) (54) 
The trajectory starts at the metastable well x{itq) = with the boundary conditions 



\ / ITQ 

For a symmetric pulse £(—t) = £(t) the velocity dx/dt = at t = 0, at this point the 



particle escape from under the barrier, and this terminates the integration in Eq. 53. The 



parameter To has a meaning of under-barrier traversal time |18| and can be found from the 
equation 

/•TQ 

£ T + / dr£{iT) = ^2m{V - E) (56) 
Jo 



Eq. [56| is equivalent to the condition p = following from Eq. |TT]. Since the time is imaginary 
the function A can be called the Eucledian action 



A = 2(V - E)r - -^r 3 / rdr \ dr^in) / dr / dr^in) 



3m m j j in 

The outgoing particle has the energy E + 5E, where 



o 



(57) 



5E = V _ E _ (£ Q + £( )) + I d r ^ d7i£(iTi)) (58) 

After escape the barrier an action of the non-stationary field on the particle can be omitted 
since it is determined by the parameter £/£q which is much smaller than one governing 
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the particle under the barrier and defined by the conditions (|38|). In the absence of a non- 
stationary pulse the energy of outgoing particles has the same value (5E = 0) and the 
Eucledian action equals Aq determined by the conventional WKB formula 

A (E) = ^(V-E)t 00 (59) 

where Too is given by Eq. [17| and has a meaning of the under-barrier traversal time in the 
stationary case. The condition 9 = tqq, which is a particular case of Eq. |3], sets some 
threshold energy 

E T = V-^ (60) 
2m 



As follows from Eqs. [57] and |58], the intensity and the energy of outgoing particles strongly 
depends on whether the initial energy E bigger (r 00 < 6) or smaller {6 < r 00 ) than E T . At 
Et < E the effect of the non-stationary signal on tunneling is weak and increases only in 
the vicinity of Et 



A = A 1 - 



38 1 \ 5E 28 1 



(n - 1)2^ (1 _ T 00 /6) n - 2 J ' V-E (n- 1)2"£ (1 - ^/ef' 1 

(61) 

At low energies E < Et the situation is very non-perturbative 

A = A (E T ) + 2(E T - E)9; 5E = E T — E (62) 

what coincides with the exponent in Eq. |46| obtained by a direct solution of the Hamilton- 
Jacobi equation. The energy dependence of the Eucledian action is shown in Fig. 4. This 
type of scenario of barrier penetration is shown schematically in Fig. 2 in the case of particle 
flux on the barrier. 



X. SEPARATION OF QUANTA ABSORPTION AND TUNNELING 



For a monochromatic field 8 U of frequency ui the total probability of penetration through 
a barrier can be approximately written as a product of two probabilities: absorption of iV 

16 



quanta and tunneling (see the comment to Eq. 1) 

exp (-A Q (E + uN)) = exp {-A{u, N)) (63) 



dt \£ ojt 00 J 

For the pulse (1) the amplitude £ UJ should be substituted by £{u6) n ~ 1 exp(— u6) and the 
effective action becomes of the form 

Afr,, JO = ME + -,A0 + 2 in ( £M) „. 2 t p( _^ ) (64) 

Semiclassical approximation corresponds to some optimum choice of uj in a continuous spec- 
trum of the pulse and the number of quanta of this optimum frequency N, which provide a 
minimum of A(uj, N). The condition dA(uj, N)/dN = gives the following relations 

8E = uN={E T -E)- 2B(y ~ E K ( So 



ur m ~ \£(lu9>- 2 
2{V-E)(, 9*\ u f £ ;i 

-2 

00 - 



A = w + w . el( _ Vl (65) 



The further minimization, with respect to uj, give an infinite (in this approach) value of uj 
indicating the logarithmic terms in Eqs. 65 to be small and hence Eqs. 65 coincide with the 
result d|). 

Let us consider another example, when such simple approach also give a correct (with 
an exponential accuracy) decay rate. Suppose the Gaussian pulse 

£{t) =£exp{-Q 2 t 2 ) (66) 

acts on a particle in the stable potential well (f§) with £q = 0. Than in the effective action 



A^N)=2N\^ m y- E) (67) 

one should put £ UJ —>■ £ exp(— oj 2 /Q 2 ), according to the Fourier harmonic of the pulse, and 
N = (V — E)/u, since in this case there is no tunneling and a particle should reach the top 
of the barrier. This leads to the relation 
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The minimization of this expression with respect to u gives the optimum value of A 

a = 2(v-e) Agv^EESV 72 (69) 

the optimum pulse frequency, and the optimum number of absorbed quanta 

■)Oli toy/miy - E) \ 1/2 M V-E 
u op t = 2S2 In — - ; N opt = (70) 

C I ^opt 



The result ([69]) coincides with calculation of the Eucledian action for the pulse (|66| ) by 
semiclassical methods developed above (calculations are not put in this paper). One can see 
from here, the decay rate under action of a non-stationary pulse can be calculated with an 
exponential accuracy on the base of simple arguments of optimum frequency and number 
of quanta. This approach of separation of quanta absorption and subsequent tunneling, 
described in this Section, works only for a potential V(x) which is not an analytical function 
of the variable x like the potential (Q). In this case one can use an interpretation of quanta 
absorption at some point x (position of singularity of V(x) on the real axis). When V(x) is 
an analytical function there is no such particular point, the situation is more complicated, 
and the method of simple separation of absorption and tunneling does not work, since the 
quantum interference of these processes becomes very non-trivial. The case of analytical 
potential is considered in the next Section. 

XI. WEAK NON-STATIONARY SIGNAL 

Let us go back to an analytical potential barrier V(x). When £ (t) = the contour C is 
reduced to the contour Cq shown in Fig. 3, which consists of the vertical part between the 
points ±z Im i and the horizontal semi- infinite lines at Ret < 0. In this static case 

Imt = ,/f / - dX (71) 
V 2 J y/V(x)-E 

where the integration goes between two classical turning points determined by the relation 
V(x) = E. Eq. [51] at £(t) = determines the conventional WKB exponent by means of the 



unperturbed Lagrangian L c 



A = -i J dtL = 2v / 2m J dx^/V{x) - E 



m ( dx 



2 



L °=2 \W) - V M + E (72) 

Here xo(t) is the classical trajectory determined at all t by Eq. |)2|. The small pulse £ (t) 
results in the perturbed trajectory Xo{t) + 5x(t). The perturbation in Eq. |5l] has the form 

A = -ij c L + ^ Xo (t)£(t)+m(^6x^ (t* Q ) - m (^Sx) (*o) (73) 

The velocities dx/dt at t$ and to are real and 

l m 5x(t ) = -(^5t )(t ) (74) 



v dt 

Here 5fo is a variation of to due to the pulse given by Eq. [fl|. One can easily see that 



J c dtL = J^dtL -2i(^j (t )lm5t (75) 



Collecting Eqs. uq-75, one can obtain 



A = A + 5A ; 6A = -i f dt £{t)x {t + At) (76) 

Jc 

Here we keep the argument shift At of the unperturbed solution, satisfying Eq. [52], deter- 
mined in the way the classical turning point t\ = —At. The method of classical trajectories 
produces a minimum value of A, this means the shift At to be found from the minimization 
condition 

85A 



dAt 



(77) 



A meaning of the minimization condition (]77 ) can be clarified in the following way. Accord- 



ing to classical mechanics, the variation of the particle energy is dE/dt = £(t)dxo/dt. At 
t = to and at t = tp, when the non-stationary field is zero, energy should have the same 
values, that is 

^ o( t + At)^ o 
c W dt 
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Eq. |78] coincides with Eq. |77|. If the minimization condition (|77|) violates the particle en- 
ergy E would acquire an imaginary part. In summary of this Section, in case of small 
non-stationary field £(t) one can use the perturbation approach ([76]) with the further mini- 
mization (fT8|). 



XII. ANALYTICAL PROPERTIES OF TRAJECTORIES 



We consider the potential barrier 



V & = J / (79) 
cosh x/a 



The classical unperturbed trajectory satisfies the relations JT5|Jl6|] 



dx (t + At) _ aujsmhuj(t + At) _ ^ 2 _ 2E ^ 

dt ^/cosh 2 u(t + At) + E/(V - E) ' mo? 

and is an analytical function of the complex variable t having the branch points at t = t s , £* 
where 



i Vv + Ve . n 

t a = it s In At; r s = — (81) 

u \JV — E 2u 



Close to the branch point t s the trajectory has the form 



x (t + At) = -™ + a\Jlu)(t s - t) ^JVjE (82) 

The cut is shown in Fig. 3 by the dashed horizontal line. Now the integral, defining 8 A in 
Eq. |7j| can be calculated on base of analytical properties. 



XIII. TUNNELING PROBABILITY 

Let us choose the non-stationary pulse in the form (Q) with n = 2 and the potential barrier 
(ff5|)- Then the integrand in Eq. |76| has singularities of two types in the complex t-plane: 



t s comes from the analytical function xo(t + At) and id comes from the analytical function 
S(t). There are different positions of the contour C with respect to those singular points 
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giving rise to different branches of the extreme action. The branch, giving the minimum 
value of A, results from the position shown in Fig. 3, when the contour C goes between the 
two singularities. This very branch determines the effect. We consider here only the case 
lmt s < 9. The contour C can be deformed up and in the limit (9 — t s ) <C 9 the pole circle 
at t — id gives the main contribution 



According to Eq. [81], t s depends on At and the minimization (|77| ) produces 

M = -t^L (84) 

The moment t\ = —At is the delay time of outgoing particles from under the barrier. The 
correction 5 A in Eq. [7B| has the form 

SA= -i £a ^( 3 zY\r^ (85) 



4 s \E J V O-Ts 
Finally, Eqs. |50|, 72, [76|, and ^5| determine the maximum outgoing flux of particles tunneling 
through the potential barrier ( |79D in the case of a weak non- stationary signal (§) with n = 2. 
Eq. [85] is analogous to the formula ([Jl]) for the triangular potential. When 9 — > Im t s the 
perturbation theory breaks down and the result becomes to be very non-linear function 
of the non-stationary pulse like Eq. ^2[ We see that the relation (0) plays a crucial role 
in physics of tunneling under non-stationary conditions. For the potential barrier (]79|) t s 
appears in a "natural" way as a result of analytical properties, whereas for a non-analytical 
potential it is determined by the time of motion between a turning point and a point of the 
non-analytical singularity of the potential (too in the case of the triangular potential (|j)). 

XIV. DISCUSSIONS AND CONCLUSIONS 

In the case of the triangular potential barrier the semiclassical theory in a non- stationary 
case is constructed based on classical trajectories in the complex time. This becomes possible 
since for such barrier non-semiclassical corrections can be calculated exactly. The conditions 
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of applicability of the semiclassical approach (softly varying and not very small pulse) are not 
surprising. Much less trivial matter is the way of construction of the trajectory method. The 
condition (|3]) of coincidence of singularity of a trajectory and of a field leads to unexpected 
conclusions. The threshold energy Et, set by the condition (|3]), divides the all incident flux 
of particles by two groups: (i) particles with E > Et passes the barrier practically like in a 
stationary case, (ii) particles with lower energies E < Et are affected strongly by the pulse 
which collects all particles after passing the barrier at the same energy Et regardless their 
initial energy. This enables to use short pulses, adapted to the necessary energy level by 
the relation (|3|), to get outgoing particles collected at this level at some moment of time. 
In other words, adapted signals can be used for a selective control of tunneling: in solids 
and molecules it excites only some particular bonds leaving other bonds non-excited. There 
is an obvious advantage of use soft signals for control of tunneling. The hard pulse of the 
type (0) cannot control tunneling selectively, since it kicks up all particles to the top of 
the barrier. When the pulse (0) serves as an envelope for the monochromatic signal of 
the frequency Q ~ 1/V it provides in a non-selective manner equal fi-shifts of all escaping 
particles with respect to their incident energies. It is remarkable, the soft Gaussian signal 
£ (t) = £exp(— 1 2 /6 2 ) cannot produce a strong enhancement of tunneling when its amplitude 
is much less than a static field of a potential barrier. This is an indication of importance of 
the analytical structure of the signal This can be understood in another way looking 
at Eq. 1, where the Fourier harmonics of the Gaussian signal £q ~ £ exp(—Q 2 8 2 /4) do not 
result in divergence in fl. 

The static electric field in solids and molecules can be estimated as £q ~ 10 7 V/cm, a 
typical pulse width is in the range of tens of femtoseconds, and the amplitude of the electric 
field of the pulse can be chosen as £ ~ 10 4 — 10 5 V/cm, what is reachable in experiments. 

The developed theory is also applicable to quantum mechanical over-barrier reflection 
of particles. The reflection of classical waves (electromagnetic, hydrodynamic, etc.) from a 
spatially-smooth medium also may be described by the above theory, when the medium is 
influenced by an adapted signal. 
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FIGURES 




FIG. 1. A plot of the imaginary part of the action for the case of the triangular potential barrier 
at t = 0. The classical position of a particle before tunneling is x = and after tunneling is x = x\ 
(a classical turning point). In absence of a non-stationary pulse the two branches go over into 
conventional increasing and decreasing WKB branches. 
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FIG. 2. A snap-shot of a particle motion through the potential barrier under the action of a 
non-stationary pulse at a moment of the maximum tunneling flux. E max is the maximum (in time) 
energy of escaped particles. All particles with the incident energies E < Et are collected at the 
threshold level Et after passing the barrier. A motion of particles with E > Et is violated a little. 
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FIG. 3. The contours of integration C and Co (in absence of a pulse) are shown in the plane 
of complex time. iO is the position of singularity of the non-stationary pulse and t s is the branch 
point singularity of the classical trajectory. The cut is denoted by the dashed horizontal line. 
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FIG. 4. The energy dependence of the exponent A{E) (solid line) which determines the max- 
imum (in time) of the tunneling probability W max ~ exp(— A(E)). The dashed line is the plot of 
the exponent Aq{E) in the absence of a pulse. It merges the solid line at E > Et- 
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